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PACS 73 . 63 . -b - Electronic transport in nanoscale materials and structures 

Abstract. - We report a theoretical low- field magnetotransport study unveiling the effect of pseu- 
dospin in realistic models of weakly disordered graphene-based materials. Using an efficient Kubo 
computational method, and simulating the effect of charges trapped in the oxide, different magne- 
toconductance fingerprints are numerically obtained in system sizes as large as 0.3/xm 2 , containing 
tens of millions of carbon atoms. In two-dimensional graphene, a strong valley mixing is found to 
irreparably yield a positive magnetoconductance (weak localization), whereas crossovers from pos- 
itive to a negative magnetoconductance (weak antilocalization) are obtained by reducing disorder 
strength down to the ballistic limit. In sharp contrast, graphene nanoribbons with lateral size as 
large as lOnm show no sign of weak antilocalization, even for very small disorder strength. Our 
results rationalize the emergence of a complex phase diagram of magnetoconductance fingerprints, 
shedding some new light on the microscopical origin of pseudospin effects. 



In 2004, the report on graphene discovery pQ has sparked a great scientific excitement 
because of both novel type of electronic excitations (so-called massless Dirac Fermions) [2] 
and a promising future of graphene-based technologies [3]. In two-dimensional graphene, 
the very peculiar nature of low-energy electronic states (encompassing a new pseudospin 
degree of freedom) yields a wealth of anomalous transport features such as Klein tunnel- 
ing [3], weak antilocalization [51|S], unconventional quantum Hall effect [7J[5], or new ways 
(supercollimation) to guide charge flows [9] . 

Besides, the nature of disorder, its effect on electronic properties, and the envisionned 
defect engineering for novel graphene electronics are currently subjects of great concern [TUj 
111] . However to date, the precise relationship between the underlying disorder features 
and the onset of graphene unique transport properties remains debated. The question also 
applies to graphene nanoribbons which exhibit different electronic properties induced by 
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lateral size confinement effects and edge symmetries [I2l - [T6] . 

One intriguing question is the precise role of pseudospin effects in the modulations of the 
conductance fingerprints in disordered graphene based materials (both two-dimensional and 
nanoribbons). Indeed, in ordinary low-dimensional disordered metals, quantum interfer- 
ences between multiple scattering paths produce a quantum correction to the semiclassical 
(Drude) conductivity which can be either reduced (weak localization) [17] or enhanced (weak 
antilocalization-WAL) |18j . depending on the strength of spin-orbit coupling. In graphene, 
it was suggested that, even in absence of spin-orbit coupling, WAL develops if the underlying 
long range disorder preserves the pseudospin symmetry of the honeycomb geometry 016]. 
A perturbation (diagrammatic) theory was proposed, but to determine the quantum correc- 
tion of conductivity, several phenomenological parameters were introduced (intravalley and 
intervalley elastic scattering times) , out of reach from analytical considerations 16. . This phe- 
nomenological approach provides a fitting equation for comparison with experimental data, 
but its connection to the specific microscopic disorder at the origin of multiple scattering 
remains totally elusive. 

A complex phase diagram of magnetoresistance phenomena in graphene has been re- 
cently experimentally revealed, exhibiting multiple crossovers from weak antilocalization 
(negative magnetoconductance) to weak localization (positive magnetoconductance) when 
tuning the energy of charge carriers or sample temperature |19j . The understanding of the 
microscopical origin of such a complex phase diagram and its relation to a possible metal- 
insulator transition remains however obscure. Its comprehension demands to clarify the 
connection between the precise disorder characteristics and the magnetic-field dependent 
quantum interference effects (QIE) driving the conductance modulations, as well as the 
true contribution of the pseudospin degree of freedom, recently questioned by Winkler and 
Ziilicke |20j . Additionally, whereas weak localization in 2D is precursor of a true insulating 
state in the zero temperature limit, WAL points towards a robust metallic state, a fact in 
marked contrast with the usual scaling theory of localization for two-dimensional disordered 
systems [17] . 

Here, by means of an efficient quantum-transport simulation method, the weak-field 
magnetoconductance of disordered graphene-based materials is explored for the case of long- 
range disorder potential, which mimicks the effect of charges trapped in the underneath 
oxide layers. This study unveils the true contribution of pseudospin effects in realistic 
disordered graphene. In two-dimensional graphene, a transition from positive to negative 
magnetoconductance is confirmed to be related to the decay of the intervalley scattering 
contribution. Our study provides however a clarification of the microscopic origin of the 
reported experimental phase diagram [19] in terms of the dominant disorder characteristics, 
whereas weak antilocalization and pseudospin effects are found to be suppressed in the 
presence of graphene edges (graphene nanoribbons). 

We use the it electrons model of graphene well described by the tight-binding Hamiltonian 
T-L = J2 a V a \a)(a\ + 70 12t a ,p) e~ %VaP with 70 = -2.7 cV the hopping integral between 

the nearest neighbor carbon atoms with distance a — 0.142 nm, whereas V a are the onsitc 
energies of 7r-orbitals and include the (long range) disorder potential profile. The effect of 
the magnetic field is introduced through the Peierls phase [24], with a flux per hexagon 
being equal to <ft = § A • d\ = h/e ^ hexagon fap- We implement a convenient gauge where 
Shoxagon fafi can take integer multiples of l/(N x N y ) to vary the magnetic field strength 
and where N x and N y = N x + 1 give the sample size. The total density of states in clean 
graphene at zero field and at B = 10 Tesla is found to exhibit the well-known \JnB Landau 
level structure in perfect agreement with theory [25| . We use a common model to describe the 
long-range disorder |21II22) , which is given by the contributions from Nj impurities randomly 
distributed at (among N sites of the system) with V a = Yli=i £ i exp(— |r Q — Ti\ 2 / (2£ 2 )), 
where £ defines the effective range, and e, are chosen at random within [— W/2, W/2] (70- 
unit), such that W gives the strength of the local potential profile (for an illustration seed] 
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Fig. 1: (a) Time-dependent diffusion coefficient at Dirac point for several impurity densities and 
W = 2.(b) Density of states of pristine graphene at zero and at 10 Tesla. (c) Schematic of the 
disordered potential in graphene due to long-range Coulomb impurities. 

(c)). Different random configurations of graphene samples with same size, £, W, and rij 
Y, Y constitute an ensemble with given disorder strength. Following prior studies [2T)l22j. 
we fix £ = 3a = 0.426 nm as a typical value for a long-range potential, but vary W to 
describe different screening situations. 

The (magneto)-transport properties of large and disordered graphene systems are sim- 
ulated by using an efficient order-N Kubo method [23]. From the study of the quantum 
wave-packet dynamics, the elastic mean free path £ e (E) can be first derived from the time 
dependence of the diffusion coefficient D(E,t) 

D(E t) = 1 EnjjWg - H)(W(t)xU(t) - x)2\* n ) 

1 E m (*™\5(E-H)\*rn) 

with x the position operator along the transport direction (at), in the Schrodinger rep- 
resentation. The behavior of D(t) is related to the time evolution of \^ n ) described by 
the operator U(t) = II n lj exp(z7JAi/ft) with At the chosen time step. An efficient Cheby- 
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shev polynomial expansion method is used for U (t) |23| . The calculations are performed 
for several initial wavepackets and for total elapsed computational times t — 9ps, with 
At = 2.44fs. The system size is about (0.5 x 0.6)/im 2 (with about 1.15 x 10 7 carbon atoms) 
and periodic boundary conditions are applied. The Kubo conductivity is finally computed 
as a(E,t = N t At) = e 2 p(E)D(E,t)/4, with p(E) the total density of states and the mag- 
netoconductance is calculated for the total elapsed time. 

In the case of disordered graphene nanoribbons (GNRs), we adopt the Landauer-Biittikcr 
approach (within the Green functions formalism) to evaluate the transmission coefficient 
T(E) and the related conductance. By performing a scaling analysis of these quantities [15] 
(averaged over a large number of different disorder configurations), we can extract the elastic 
mean free path £ e . 

Q] (a) shows the evolution of the wavepacket dynamics (through time-dependent diffusion 
coefficients) at the Dirac point (E=0) for W = 2 and increasing impurity density (m — 
0.125%, 0.25%, 0.5%). The diffusion coefficients are found to reach a saturation regime after 
about lps, and then to remain almost constant up to very long times, indicating a diffusive 
regime. At long elapsed propagation times, the saturation of the diffusion coefficient is 
followed by a time-dependent decay (for all energies and considered impurity densities), 
which pinpoints an increasing contribution of QIE and weak localization phenomena. 

The maximum value of D(t) — Z? m ax allows for the evaluation of the elastic mean free 
path £ e = D Wl3X /2vf (cf. H]) where vp — 8.7 x 10 5 ms _1 is the Fermi velocity. It is seen that 
t e ~ 1 /rii in agreement with Fermi's Golden rule. The dependence of the elastic scattering 
time r e = £ e /vF on the charge density n follows r e ~ f° r sufficiently small n and W 

(not shown here). 

The effect of an applied magnetic field is shown in [3] Top and bottom panels give the 
time-dependent diffusion coefficients at varying magnetic field and the magnetoconductance 
Aa(B) = a(B) — <r(B = 0), respectively. For W = 2 {3£a)), a steady suppression of QIE- 
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Fig. 3: Top panels: diffusion coefficient at the Dirac point for various magnetic fields and W, 
rii — 0.125%, £ = 0.426 nm. Bottom panels: Aa(B) for four different Fermi level positions 
{E-dp = 0,Ei = 0.049 eV, En = 0.097 eV and Em = 0.146 eV). Dashed lines are fits as explained 
in the text. For W = 1.0 in (f) data and fit at -Bdp have been rescaled by 0.3 for clarity. 



related conductance corrections is clearly seen with increasing B up to B = 0.274: T at which 
the decay of D(t) has been fully suppressed. This is in accordance with the weak-localization 
phenomenon, as confirmed by the positive magnetoconductance reported in panel (d) for 
four selected energies (averages over 32 different configurations are performed). In contrast, 
the case W — 1.5 shows a clearly different behavior as seen in [3Jb) and (e) respectively 
for D(t) and Aa(B) (for which 64 configurations have been averaged). Indeed, a change 
in the sign of Aa(B) is observed, in agreement with a crossover from weak localization to 
weak antilocalization, which depends not only on the field strength, but also the considered 
energy of charge carriers. We note that an estimation of the elastic mean free path gives 
£ e ~ 9 — 20nm, which is much smaller than our sample size. 

For the case W = 1, we first observe that at zero magnetic field, D(t) does not reach a 
saturation regime within our total computational time (dashed black curve in[3](c)). With 
increasing the magnetic field, D(t) is found to be reduced and no onset of weak localization 
manifests. The corresponding Aa(B) curves displayed in panel (f) all exhibit a negative 
magnetoconductance, whose connection to a WAL regime, as discussed here below, is how- 
ever not straightforward. 

To rationalize all those results, we first note that the numerical study reported in |22j 
suggests that as long as W > 1, the long-range potential (with £ = 0.426nm) introduces 
strong valley mixing. This is consistent with our computed positive magnetoconductance 
for the case W = 2 (0Jd)). By further decreasing the disorder strength (from W = 2 
to W — 1.5), WAL occurs as seen in [3fe) which corresponds to the decay of intervalley 
scattering contribution (see also Figure 4 of [22]). Our study thus exhibits a transition 
from weak localization (for W = 2) to WAL (W = 1.5), supporting that the enhancement 
of intervalley scattering contribution yields predominance of the weak localization effect. 
The analysis can be deepened by fitting the simulations using the phenomenological law 
Aa(B)=e 2 /irh{F(Ts 1 /T- 1 )-3F(TB 1 /(T- l +T- 1 ))}, where F{z) = lnz + ^(1/2 + z" 1 ), 
ip(x) is the digamma function and Tg 1 = 4eDB/h [6]. Here, for simplicity we restrict the 
model to a single additional elastic scattering time r* (phenomenological parameter) that 
contains both contributions of intravalley and intervalley scattering [B]. The least-squares 
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Fig. 4: (a) Inset: Transmission versus energy for a 10 nm-wide armchair ribbon (84-aGNR) with 
W = 1/2, rii = 0.1% and for different lengths (up to 1.5/xm) of the disordered region. Main frame: 
Energy-dependent mean free path for m = 0.1% (black curve) and m = 0.2% (red curve), (b) 
Schematic of a short segment of a disordered nanoribbon with long-range scatters. Bright (dark) 
spots correspond to positive (negative) values of V. 

fits (dashed lines) are superimposed to the simulated Aa(B) (symbols), taking t v = 9ps 
(our maximum computed time). For W = 2 and W = 1.5, t» ranges within [1.1 — 2.3]ps 
and [1.5 — 6.3]ps, respectively (increasing values with increasing energy), thus confirming 
the weak localization regime for lowest B, which is fully consistent with [5] (r, < r„). The 
case W=l is more complicated, with a fit (with t* > t v ) only possible at lower magnetic 
fields and higher energies (En and Em) and out of reach for energies nearby the Dirac 
point. Despite the obtained negative magnetoconductance, the interpretation in terms of 
weak antilocalization is actually not obvious. Indeed, as seen on the diffusion coefficient 
behavior at zero magnetic field, the saturation regime is not reached within our maximum 
computed time, although the conduction clearly departs from a ballistic regime. By using 
£ e = D max /2vF, with an extrapolated maximum diffusion coefficient of > 2500cto 2 s~ 1 , one 
gets £ e > 330nm, which is in the order of our graphene longitudinal size. Accordingly, the 
transport regime is closer to a quasiballistic case, which jeopardizes the interpretation in 
terms of WAL. 

To analyze the contribution of sample geometries (edges), we finally simulate numerical 
transport in graphene nanoribbons [T2J[T5]. As a typical case, we focus on a 10 nm-wide 
armchair ribbon (84-aGNR) with impurity densities »j = 0.1%, 0.2%, £ = 0.426 nm and 
W = 1/2. The results for T(E) are averaged over 100 disorder configurations (JU-inset). At 
the low fields (up to 0.5 T) considered here, the bandstructure does not change significantly 
with respect to the zero field case. Accordingly, the impact of the field on transport is only 
determined by the modifications of QIE. 2] (main frame) shows the elastic mean free path of 
disordered ribbons as a function of electron energy. As in the case of Anderson-type disorder 
or edge defects |15j . £ e exhibits important modulations depending on the electron energy, 
with systematic decay around the energies corresponding to the onset of new subbands. At 
low energies, £ e can even reach a few microns for the considered disorder values. 

E shows the evolution of the quantity AG(B) = G(B) - G(B = 0) (for m = 0.1%) for 
the same set of ribbon lengths considered in|4j and at two different selected energies (dashed 
vertical lines of HJ). Note that the small value of W = 0.5 considered here would prohibit 
valley mixing in the two-dimensional case. The ratio L/£ e gives a suitable criterion to classify 
the transport regime as quasi-ballistic (when L/£ e < 1) and diffusive or localized (when 
L/£ e > 1). When L/£ e > 1, a markedly positive magnetoconductance dominates (AG(B) > 
0) in agreement with the standard weak localization regime ([5]). The magnetoconductance 
variation increases for larger lengths, giving a stronger field-driven suppression of QIE. The 
results in the case rii = 0.2% are analogous. Our numerical results evidence the absence of 
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Fig. 5: Magnetoconductance for the disordered ribbon (W = 0.5) at two different selected energies 
(see 2] vertical dashed lines) and for different ribbon lengths L. Different set of curves are identified 
depending on the ratio L/£ e . 

weak antilocalization for disordered GNRs because of the strong reintroduction of intervalley 
scattering in presence of ribbon edge geometries. One notes that the application of much 
larger magnetic fields has been recently theoretically predicted and experimentally found to 
yield very large magnetoresistance signals. [2fjll2T] 

In conclusion, in two-dimensional disordered graphene the crossover between weak local- 
ization and weak antilocalization was theoretically demonstrated to be driven by the nature 
of long-range Coulomb scattering potential. In contrast, predominance of weak localization 
effects was obtained in graphene nanoribbons. These results allow to bridge the nature of 
disorder in graphene with the localization phase diagram explored experimentally |19j and 
also open new directions for understanding pseudospin effects in graphene-based materials. 
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